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Abstract 

It is important to provide guidance on whether CP violation may be measurable 
in top-quark production at the Large Hadron Collider. The present work extends an 
earlier analysis of the non-supersymmetric Two-Higgs-Doublet Model in this respect, 
^ . by allowing a more general potential. Also, a more comprehensive study of theoretical 

and experimental constraints on the model is presented. Vacuum stability, unitarity, 
direct searches and electroweak precision measurements severely constrain the model. 
We explore, at low tan /3, the allowed regions in the multidimensional parameter space 
that give a viable physical model. This exploration is focused on the parameter space 
of the neutral sector rotation matrix, which is closely related to the Yukawa couplings 
of interest. In most of the remaining allowed regions, the model violates CP. We 
present a quantitative discussion of a particular CP-violating observable. This would 
be measurable in semileptonically decaying top and antitop quarks produced at the 
LHC, provided the number of available events is of the order of a million. 
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1 Introduction 



The Two-Higgs-Doublet Model (2HDM) is attractive as one of the simplest extensions 
of the Standard Model that admits additional CP violation [1-3]. This is an interesting 
possibility, given the unexplained baryon asymmetry of the Universe [4,5], and the pos- 
sibility of exploring relevant, new physics at the LHC [6]. In particular, the model can 
lead to CP violation in ti production, a process which has received considerable theoretical 
attention [7-9], since it will become possible to severely constrain or even measure it. 

CP violation can be induced in ti production at the one-loop level, by the exchange of 
neutral Higgs bosons which are not eigenstates under CP. This effect is only large enough 
to be of experimental interest if the neutral Higgs bosons are reasonably light, and have 
strong couplings to the top quarks. 

Within the 2HDM (II), where the top quark gets its mass from coupling to the Higgs 
field $2 [10] (see sect. 13.51) . the condition of having sizable Hti couplings forces us to 
consider small values of tan/3. A first exploration of this limit was presented in [11]. In 
that paper, the general conditions for measurability of CP violation in gg — ► ti at the 
LHC [8] were found to be satisfied in a certain region of the 2HDM parameter space. In 
addition to having small tan /3, in order to have a measurable signal with a realistic amount 
of data (of the order of a million ti events), it was found necessary that the lightest neutral 
Higgs boson be light, and that the spectrum not be approximately degenerate. In fact, it 
was found that in the most favourable observable considered, the effect would not reach 
the per mil level unless there is one and only one Higgs boson below the ti threshold, and 
that tan /3 is at most of order unity. We here extend the analysis of [11] to the more general 
case, allowing the most general quartic couplings in the potential. 

At small tan f3, also certain Yukawa couplings to charged Higgs bosons are enhanced. 
Such couplings contribute to effects that are known experimentally to very high precision. 
In particular, at low tan/3 the B®-B® oscillation data and the effective Zbb coupling, 
measured via [12, 13] severely constrain the model, whereas the b — > 37 data [14] 
constrain it at low M#± . Furthermore, the high-precision measurement of the W and Z 
masses, as expressed via p [13] constrains the splitting of the Higgs mass spectrum. Unless 
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there are cancellations, the charged Higgs boson can not be very much heavier than the 
lightest neutral one, and the lightest neutral one can not be far away from the mass scale 
of the W and the Z [15]. Also, the lightest one is constrained by the direct searches at 
LEP [16,17]. We shall here study the interplay of these constraints, and estimate the 
amount of CP violation that may be measurable at the LHC in selected favorable regions 
of the remaining parameter space. 

An important characteristic of the 2HDM (as opposed to the MSSM [18-22]) is the 
fact that, at the level of the mathematics, the masses of the neutral and the charged 
Higgs bosons are rather independent (see sect. However, the experimental precision on 
Ap (see sect. 14.31) forces the charged Higgs mass to be comparable in magnitude to the 
neutral Higgs masses. Another important difference is that whereas small values of tan /3 
are practically excluded in the MSSM [23], in the 2HDM, which has more free parameters, 
they are not. 

For a recent comprehensive discussion of the experimental constraints on the 2HDM 
(though mostly restricted to the CP-conserving limit), see [24] and [25]. The latter study, 
which considers the CP-conserving limit, concludes that the model is practically excluded, 
with the muon anomalous moment being very constraining. However, the interpretation 
of the data is now considered less firm, and furthermore, that study focuses on large tan f3, 
and is thus less relevant for the present work. 

We present in sect. [2] an overview of the 2HDM, with focus on the approach of ref. [11], 
and outline the present extensions. In sect. [3] we discuss the model in more detail, in 
particular the implications of stability and unitarity, and review the conditions for having 
CP violation. In sect. H] we discuss various experimental constraints on the model, with 
particular attention to small values of tan/3. In sect. E]we present an overview of allowed 
parameter regions, also restricted to small tan/5. In sect. |6]we discuss the implications of 
the model for a particular CP-violating observable involving the energies of positrons and 
electrons from the decays of t and i produced in gluon-gluon collisions at the LHC. Sect. [7] 
contains a summary and conclusions. 
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2 Review of the Two-Higgs-Doublet Model 



The 2HDM may be seen as an unconstrained version of the Higgs sector of the MSSM. 
While at tree level the latter can be parametrized in terms of only two parameters, con- 
ventionally taken to be tan/3 and Ma, the 2HDM has much more freedom. In particular, 
the neutral and charged Higgs masses are rather independent. 

Traditionally, the 2HDM is defined in terms of the potential. The parameters of the 
potential (quartic and quadratic couplings) determine the masses of the neutral and the 
charged Higgs bosons. Alternatively, and this is the approach followed here and in ref. [11], 
one can take masses and mixing angles as input, and determine parameters of the potential 
as derived quantities. This approach highlights the fact that the neutral and charged sectors 
are rather independent, as well as masses being physically more accessible than quartic 
couplings. However, some choices of input will lead to physically acceptable potentials, 
others will not. This way, the two sectors remain correlated. 

In addition, the 2HDM neutral sector may or may not lead to CP violation, depending 
on the choice of potential. We shall here consider the so-called Model II, where w-type 
quarks acquire masses from a Yukawa coupling to one Higgs doublet, $ 2 , whereas the 
rf-type quarks couple to the other, $x. This structure is the same as in the MSSM. 



The amount of CP violation that can be measured in ti production was related to the Higgs 
mass spectrum and other model parameters in [11]. In that paper, the Higgs potential 
studied was parametrized as [26] 



2.1 The approach of ref. [11] 



V = -^(^l^x) 2 + ^(«^t^ 2) 2 + a 3 (^^ 1 )(^>^ 2 ) + A 4 (^ t 1 ^ 2 )(^^ 1 ) 
+i [A 5 ($i$2) 2 + h.c. 

-- {m^SlSi) + [m? 2 ($l$ 2 ) + h.c] + T74($t$ 2 ) j _ 



(2.1) 



Expanding the Higgs-doublet fields as 




(2.2) 
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and choosing phases of such that v\ and v 2 are both real [27], it is convenient to define 
773 = — sin Pxi + cos @X2 orthogonal to the neutral Goldstone boson G° = cos ftxi + sin 0X2- 
In the basis (771, 772, 773), the resulting mass-squared matrix M 2 of the neutral sector, can 
then be diagonalized to physical states (Hi, H2, H$) with masses Mi < M2 < M3, via a 
rotation matrix R: 









H 2 ) 




(2.3) 


hJ 


- R : ■ 





satisfying 



and parametrized a; 



i?.M 2 i? T = M 2 = diag(M 1 2 , M 2 , M 3 2 ), (2.4) 



1 \ / cos 0:2 sina^X / cosai sinai 

R = R 3 R 2 Ri = I cos as sin II 1 II — sin a\ cos a\ 
^0 —sin 0:3 cos as J \— sin a 2 cos 0^2/ \ 1 

ClC 2 SiC 2 s 2 \ 

(ci s 2 s 3 + si c 3 ) C1C3-S1S2S3 c 2 s 3 (2.5) 
-C1S2C3 + S1S3 -(cis 3 + sis 2 c 3 ) c 2 c 3 / 

with q = cosctj, Sj = sinctj. The rotation angle ai is chosen such that in the limit of no 
CP violation (s 2 — > 0, S3 — > 0) then a>i — > a + -^k, where a is the familiar mixing angle 
of the CP-even sector, and the additional ^tt provides the mapping Hi <-> h, instead of H 
being in the (1, 1) position of -M 2 iag , as is used in the MSSM [10]. 

While the signs of rji and x% are fixed by our choice of taking the vacuum expectation 
values real and positive [27], the phase of Hi has no physical consequence. One may 
therefore freely change the sign of one or more rows, e.g., let Ru —Rn (see sect. 13.1.11) . 

Rather than describing the phenomenology in terms of the parameters of the potential 
(12. ip . in [11] the physical mass of the charged Higgs boson, as well as those of the two 
lightest neutral ones, were taken as input, together with the rotation matrix R. Thus, the 
input can be summarized as 

Parameters: tan/3, (Mi, M 2 ), (M H ±, fi 2 ), (ai, a 2 , a 3 ), (2.6) 

where tan/3 = v 2 /vi and /j, 2 = v 2 v, with v = Re m\ 2 / (2viv 2 ) and v = 246 GeV. 



1 In ref. [11], these angles were referred to as (a, at,, a c ) «-> (aj., a%, a^). 
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This approach provides better control of the physical content of the model. In partic- 
ular, the elements R13 and R23 of the rotation matrix must be non-zero in order to yield 
CP violation. For consistency, this requires ImAs and l~mm\ 2 ( as derived quantities) to be 



non-zero. 



2.2 The general potential 

For the potential, in this study, we take 

V = Expression flUD + { + A r( $ 2 $ 2)] ($1* 2 ) + h.c.} . (2.7) 

The new terms proportional to Xq and A7 have to be carefully constrained, since this 
potential does not satisfy natural flavour conservation [28], even if each doublet is coupled 
only to up-type or only to down-type flavours. 

The various coupling constants in the potential will of course depend on the choice of 
basis ($1, $2)- Recently, there has been some focus [29] on the importance of formulating 
physical observables in a basis-independent manner. Here, we shall adopt the so-called 
Model II [10] for the Yukawa couplings. This will uniquely identify the basis in the ($1, $2) 
space. 

Minimizing the potential ( \2.7\\ . we can rewrite it (modulo a constant) as 



V 



Ai 
2 



+ 



A, 



.21 2 



+ A 3 ($ t 1 $l)($ 2 $2) + A 4 ($1$2)($S$1 



+ \ -A 5 ($I$ 2 ) 2 + [A 6 ($5$i) + A 7 ($S$ 2 )] ($I$ 2 ) + h.c. 



~ [Re A34567 - 2v] [vf ($t$i) + u?($J$ 2 )] - Vl v 2 Re [A 6 ($ t 1 $ 1 ) + A 7 ($^ 2 )] 



—viv 2 [2u Re ($1$ 2 ) - Im A 567 Im ($I$ 2 )] • 
Here and in the following, we adopt the abbreviations 
A345 — A3 + A4 + A5, A34567 
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A345 H A6 H A7, A567 

v 2 Vt 



X 5 + -Xe + -X 7 . (2.9) 

V 2 Vi 



The mass-squared matrix M 2 of (12.41) . corresponding to the neutral sector of the po- 
tential, is found to be 

M 2 U = v 2 [clX 1 + s 2 ^+^Re(3cf i X 6 -slX 7 )], 
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Mj 2 



v 2 [s 2 pX 2 + CpV+ ^j-Re(-c^A 6 + 3sp\ r j 



M 



2 

33 



w 2 Re [-A 5 



'7? 
1 



.M 2 



12 



M\ z 



% A 6 + . 



f 2 [c /3 s / 3(ReA 3 45 -v) + §Re(c§A 6 + s 2 X 7 )], 
-\v 2 lm [3/3X5 + 2c^A 6 ], 
-±v 2 Im [c p X b + A 7 ], 



(2.10) 



with M% = M%. 

Here, compared with the potential ( 12. ip . we have two more complex parameters, Ag and 
A7 (four new real parameters), but rather than those, we take as additional parameters 
M 3 , ImA 5 , ReA 6 and ReA 7 . Thus, the input will be 

Parameters: tan/?, (Mi, M 2 , M 3 ), (M#±, /i 2 ), (a>i, a>2, a 3 ), ImAs, (ReA6, ReA 7 ). 

(2.11) 

3 Model properties 



We want to explore regions of parameter space where there is significant CP violation. In 
order to do that, we need to map out regions in the {a±, «2, ^3} space where the model is 
consistent (figures are presented in sect. EJ). 
From eq. (12 ,4p . it follows that 



Ml = Y. Rk i M * R *r 



(3.1) 



Here, it is evident that the signs of the rows of R play no role. 
Comparing the expressions (13. ip with (I2.10p . invoking also 



M 2 H± = fx 2 - V(A 4 + Re A 



567j 



we can solve for the A's. In particular, it follows from (12.101) that 



Im A 6 = 

Im A 7 = 

S /3 



4.M 2 3 + ^ImA £ 



4^23 + yl- A, 



(3.2) 



(3.3) 
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3.1 Symmetries 

By exploiting certain symmetries of the rotation matrix R, we can reduce the ranges of 
parameters that have to be explored. 

3.1.1 Transformations of the rotation matrix 

The rotation matrix R is invariant under the following transformation; 

A : oti — > 7r + oti, a 2 ~~ > vr — ^2; ct 3 — > 7r + a 3 ; (3-4) 

which leaves its elements unchanged. 

Another class of transformations are those where two rows of R (i.e., physical Higgs 
fields) change sign, as discussed in sect. 12.11 The transformations are [11]: 

Bl : a% — > 7r + oji, a 2 — > 7r - a 2 , «3 fixed : — > i?n, i?2i — > — -Km, #3i — > — -R31) 

B2 : ai fixed, a 2 -> vr + a 2 , a 3 -> -a 3 : i? M -> --Rh, i?2i -> #2i, #3* — > -^3i, 

B3 : ai — > + ai, a 2 — > — at 2 , a 3 — > — a 3 : — > — ifcu, _R 2 i — > -#2*, #3i — ► #3i- 

(3.5) 

Actually, any one of these is a combination of the other two. For example, the transfor- 
mation B3 is the combination of Bl and B2. Other transformations exist that will yield 
the same symmetries, but they will be combinations of one of these three transformations 
followed by the transformation A. In total we have 6 different transformations that yield 
symmetries of type B. 

The third class of transformation we consider are those where two columns of R change 
sign. These transformations are: 

CI : an — > 7r - «i, a 2 — > vr + a 2 , a 3 fixed : Rji — > Rji, R j2 — > — -Rj2, -Rj3 — > --R73, 

C2 : «i — > —ai, a 2 — > 7T + a 2 , a 3 fixed : -Rji — > —Rji, Rj 2 — > -Rj 2 , -R73 — >■ —Rj3, 

C3 : «i — > 7T + «i, a 2 fixed, a 3 fixed : Rji — > —Rji, Rj 2 — * —Rj2, Rjz ~^ Rj3- 

(3.6) 

The transformation C3 is the combination of the transformations CI and C2. Other 
transformations exist that will yield the same symmetries, but they will be combinations 
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of one of these three transformations followed by the transformation A. In total we have 6 
different transformation that yield symmetries of type C. 

Under transformations of type A and B, the resulting mass-squared matrix Ai 2 = 
i? T A4^ iag -R will be invariant. We make use of this fact along with the symmetries A, Bl 
and B2 to reduce the parameter space under consideration to 

- tt/2 < { ai , a 2 , a 3 } < tt/2. (3.7) 

Under transformations of type C, the mass-squared matrix will not be invariant, some 
of its non-diagonal elements will change sign while the rest are unaltered. 

Ci: M 2 12 ^-M 2 12 , Mj 3 ^-M 2 13 , M 2 23 -*+M 2 23) 
C2: M\ 2 ^-M\ 2 , M\ 3 ^+M\ 3 , Mj 3 ^-M 2 23 , 
C3: M 2 12 ^+M 2 12 , M 2 3 ^-M 2 3 , M\ 3 ^-M\ 3 . (3.8) 

While a change of the sign of M.\ 2 implies changes in the physical content of the model, 
a change of sign of M.\ 3 and/or Ad 23 can be compensated for by adjusting the imaginary 
parts of A5, A6 and A7. Thus, the most interesting transformation among the set (13. 8ft is 
C3. 

The transformation B3 • C3 is physically equivalent to C3 since transformations of type 
B leave the mass-squared matrix invariant: 

B3-C3 : a x fixed, a 2 -> -a 2 , a 3 -»• -a 3 : M\ 2 -^ +M 2 12 , M 2 13 -^ -M 2 13 , M\ 3 -^ -M 2 23 . 

(3-9) ' 

When IH1A5 = 0, it follows from (13. 3p that a sign change of M.\ 3 and M.\ 3 can be 
compensated for by sign changes of ImAg and IrnXj. These signs play no role in the 
discussion of stability (see Appendix A) and unitarity [30]. We shall therefore, when 
discussing the case Im A5 = (sects. 15.11 and 1572]) . make use of (13.91) to restrict the angular 
range from (13.71) to the smaller 

- tt/2 < {ax, a 2 } < tt/2, < a 3 < tt/2. (3.10) 

When Im A5 7^ we need to consider the angular range as given in (13 .7\\ . 



9 



3.1.2 Inversion of tan/3 

The Higgs sector is invariant under 



tan (3 <-> cot /3, a± <-> |7r — cci, 



accompanied by 



Ai <-> A 2 , A3, A 4 <->• A 3 , A 4 , 

A 5 <-> A* A 6 <-> A*. (3.12) 

This is just the symmetry between <3>! and $ 2) and will be violated by the introduction 
of Model II Yukawa couplings, which distinguish between the two Higgs doublets, i.e., 
between tan/? and cot/3. 

3.2 CP violation 

In general, with all three rotation-matrix angles non-zero, the model will violate CP. How- 
ever, in certain limits, this is not the case. In order not to have CP violation, the mass- 
squared matrix must be block diagonal, i.e., one must require 

M 2 13 = M 2 23 = 0. (3.13) 

Thus, CP conservation requires 

Mj 3 = M 2 R U R 13 + M 2 R 21 R 23 + M 2 R 3l R 33 = 0, 

M 2 23 = MfR l2 R 13 + M\R 22 R 23 + M 2 R 32 R 33 = 0. (3.14) 

One possible solution of (13.141) is that 

Mi = M 2 = M 3 . (3.15) 

The expressions (I3.14p then vanish, by the orthogonality of R. There are additional limits 
of no CP violation, as discussed below. 
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Expressed in terms of the angles of the rotation matrix, the above elements describing 
mixing of the CP-even and CP-odd parts of M 2 take the form 



M\ z = c lC2 s 2 {Ml - s\M\ - c\Ml) + Sl c 2 c 3 s 3 (-M 2 + M 3 2 ), 

M 2 23 = cic 2 c 3 s 3 (M 2 2 - Ml) + sic 2 s 2 (M 2 - s 2 M 2 - c 2 M 3 2 ). (3.16) 

In the mass-non-degenerate case, they vanish (there is thus no CP violation) if either: 

Case I: sin2a 2 = 0, and sin2a 3 = 0, or 

Case II: cosct 2 = 0, ct 3 arbitrary. (3-17) 

Note that M\ - s§M| - c§Mf < for non-degenerate or partially degenerate masses, 
ordered such that M\ < M 2 < M 3 (where no more than two of the masses are equal). 
Thus, there are no additional CP-conserving solutions for the vanishing of this factor. The 
cases of partial degeneracy, Mi = M 2 ^ M 3 , and Mi ^ M 2 = M 3 will be discussed in 
sect. 13.71 

It is thus natural to focus on the angles a 2 and a 3 . In particular, since i?i 2 i?i 3 is 
associated with CP- violation in the Hitt coupling (see sect. 13.51) . we are interested in 
regions where | sin(2o; 2 )| is large. 

3.3 Reference parameters 

In order to search for parameters with "large" CP violation, we will assume Hi is light, 
and that M 2 is not close to Mi, as such degeneracy would cancel any CP violation. 
For illustration, as a conservative default set of parameters, we take 

Set A- J tan ' 9 = {0-5,1-0,2.0}, (M 1 ,M 2 ,M 3 ) = (100,300,500) GeV, 
6 ' \M H ± = 500 GeV, n 2 = (200 GeV) 2 , ImA 5 = 0, ReA 6 = ReA 7 = 0. l ' ' 

Here, the lightest neutral Higgs boson can be accommodated by the negative LEP searches 
[16, 17] provided it does not couple too strongly to the Z, and the charged Higgs boson 
mass is compatible with the negative LEP [13] and Fermilab searches [31] as well as with 
the B^-B^ oscillation, R b constraints [25] (see sect. 14. ip and the b — >• 57 analysis at low 
tan/5 [14]. 
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As a second set of parameters, we take 



Set B: 



tan (3 = {0.5,1.0,2.0}, (M 1 ,M 2 ) = (80,300) GeV, M 3 = {400,600} GeV 
M H ± = 300 GeV, /i 2 = {0, (200 GeV) 2 }, ImA 5 = 0, Re A 6 = Re A 7 = 0. 



(3.19) 



This set, which represents a light Higgs sector, is marginally in conflict with data (the 
combination of charged-Higgs mass and tan/3 values violate the Rb constraints by up to 
5a, see Table[2]in sect. 14. 5p . but is chosen for a more "optimistic" comparison, since it could 
give more CP violation due to a lower value of M\ (which enhances the loop integrals). 

3.4 Stability and unitarity 

A necessary condition we must impose on the model, is that the potential is positive when 
|$i| and |$2| — ► oo. This constraint, which is rather involved, is discussed in Appendix A. 
Two obvious conditions are that 



In general, the additional stability constraint is that A3 and A4 cannot be "too large and 
negative", and that [ As J , |Ae| , |A 7 | cannot be "too large". 

Furthermore, we shall impose tree-level unitarity on the Higgs-Higgs-scattering sector, 
as formulated in [30,32] (see also ref. [33]). This latter constraint is related to the per- 
turbativity constraint (A's not allowed "too large") adopted in ref. [11], but actually turns 
out to be numerically more severe. 

3.5 Yukawa couplings 

With the above notation, and adopting the so-called Model II [10] for the Yukawa couplings, 
where the down-type and up-type quarks are coupled only to $1 and $ 2 , respectively, the 
couplings can be expressed (relative to the SM coupling) as 



Ai > 



A 2 > 0. 



(3.20) 



HM : 



1 



[Rj! - inf 5 sin PR j3 ], 



cos (3 
1 



[Rj 2 - 275 cos /3Rj 3 ] = a + ia^. 



(3.21) 



HM : 



sin (3 
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Likewise, we have for the charged Higgs bosons [10] 

H + bi : — % [m 6 (l +7 5 )tan/? + m t (l - 75 )cot/3], 

2\ / 2m w 

H~tb: — % i- — [m 6 (l-7 5 )tan/3 + m t (l+7 5 )cot/3]. (3.22) 
2y2m w 

With this Yukawa structure, the model is denoted as the 2HDM (II). 
The product of the Hjti scalar and pseudoscalar couplings, 

m cos/3 
7 W =- a ~ a =-^L R R 3.23 

sin p 

plays an important role in determining the amount of CP violation in the top-quark sector. 

As was seen in ref. [11], unless the Higgs boson is resonant with the tt system, CP 
violation is largest for small Higgs masses. For a first orientation, we shall therefore 
focus on the contributions of the lightest Higgs boson, H±. (There will also be significant 
contributions from the two heavier Higgs bosons, as discussed in sect. El) For the lightest 
Higgs boson, the coupling f)3.2ip becomes 

tt - 1 r i dl 1 sin ai sin(2a 2 ) 

H\tt : sin ai cosa2 — ijs cosp sin a 2 , with X r p = o 7, 7, — > (3.24) 

sinp z tan p sin p 

where a\ and ct2 are mixing angles of the Higgs mass matrix as defined by Eqs. 02. 4ft and 



From (I3.24p . we see that low tan/3 are required for having large CP violation in the 
top-quark sector. However, according to (13.220 . for low tan/3 the charged-Higgs Yukawa 
coupling is also enhanced. Thus, for low tan /3, the Rb, AM^ [12,13] and b — * 57 constraints 
[14] force M H ± to be high. For a quantitative discussion, see sect. 14.11 

3.6 CP violation in the Yukawa sector 

We shall in sect. [6] study CP violation in the process 

pp -> ttX -> e + e~X, (3.25) 

focusing on the sub-process 

gg^ti^ e + e-X. (3.26) 
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Let the CP violating quantity of that process be given by [8, 11] 

3 

Y^R^RjsfiMj), (3.27) 
j'=i 

where f(Mj) is some function of the neutral Higgs mass Mj, in general determined by loop 
integrals. 

When the three neutral Higgs bosons are light, they will all contribute to the CP- 
violating effects. In fact, in the limit of three mass-degenerate Higgs bosons, the model 
may still be consistent in the sense that solutions can be found in some regions of parameter 
space, but the CP violation will cancel, since [cf. eq. (13.231) ] 

£tS = |^I> 2 ^3 = (3.28) 

due to the orthogonality of R. 

3.7 Degenerate limits 

The set of free parameters (12.111) permits all three neutral Higgs masses to be degenerate. 
As discussed above, in this limit there is no CP violation, by orthogonality of the rotation 
matrix R. However, in contrast to the case of Xq = A7 = studied in [11], the partial 
degeneracies are non-trivial and may lead to CP violation for certain choices of the angles 

OLi- 

Mx = M 2 = m 7^ M s . In this limit, the elements of M 2 that induce CP violation, are 

Mj 3 = R 31 R 33 (Mi - m 2 ) = c 2 c 3 ( Sl s 3 - Cl s 2 c 3 )(Ml - m 2 ), 

M\ 3 = R 32 R33(M 3 2 - m 2 ) = -c 2 c 3 ( Cl s 3 + s lS2 c 3 )(M 2 - m 2 ). (3.29) 

These both vanish, when the conditions (13.171) are satisfied, or else, when 

cosa3 = 0, ai,«2 arbitrary. (3.30) 

By orthogonality, when the two lighter Higgs bosons are degenerate, the CP violation 
(I3.27P in the top-quark sector is proportional to 

#32 R 3 3[f(M 3 ) - f(m)} ~ M 2 23 . (3.31) 
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Thus, even though the model violates CP in the limit Mi = M 2 7^ M 3 , by for example 
having M.\ 3 7^ 0, the top-quark sector would not violate CP at the one-loop level unless 

M\ 3 ~ #32^33 7^ 0. 

Mi 7^ M 2 = Ms = M. In this limit, the elements of Ai 2 that induce CP violation are 

M 2 13 = -R n Ri 3 (M 2 - Ml) = - Cl c 2 s 2 {M 2 - Ml), 

M 2 23 = -R l2 R 13 {M 2 - Ml) = - Sl c 2 s 2 {M 2 - Ml). (3.32) 

We note that these both vanish for sin(2a2) = 0, meaning a 2 = or a 2 = n/2. Thus, in 
the limit Mi 7^ M 2 = M 3 and a 2 = 0, but a 3 arbitrary, the model does not violate CP, in 
agreement with the results of [11]. 

In this limit of the two heavier Higgs bosons being degenerate, the CP violation in the 
top-quark sector is proportional tCLl 

R 12 ^[/(Mx) - f(M)) ~ M\ 3 . (3.33) 

In our parametrization, this is non-zero for 

sinai 7^ 0, sin 2^2 7^ 0, (3.34) 

but with a 3 arbitrary. 

In the more constrained model discussed in [11], the latter limits of only two masses 
being degenerate do not exist. In that case, with A6 = A7 = 0, a degeneracy of two masses 
forces the third one to have that same value. 

4 Experimental model constraints at low tan f3 

It is convenient to split the experimental constraints on the 2HDM into two categories. 
There are those involving only the charged Higgs boson, H ± , and those also involving the 
neutral ones. The former, like the non-discovery of a charged Higgs boson, the b — > 57 
constraint [14], and the B®-B® oscillations [25] do not depend on the rotation matrix R 

2 Whereas both these degenerate limits yield CP-violation in the t-quark sector proportional to .M 2 ^, 
the corresponding quantities in the 6-quark sector are proportional to M.\%- 
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and the amount of CP violation. They are given by Mh± and its coupling to quarks, 
( 13.22(1 . i.e., on tan/?. On the other hand, constraints involving the neutral ones depend on 
the details of the couplings, i.e., they depend sensitively on the rotation matrix R as well 
as on the neutral Higgs mass spectrum. We shall first review the constraints that depend 
only on the charged Higgs sector. 

In subsections I4.2H4.5I we discuss constraints on the model that depend on the neutral 
sector. For the purpose of determining these constraints, one has to generalize some pre- 
dictions for the CP-conserving case to the CP- violating case. Eqs. (14. 3ft . (14.10(1 . (14.121) . 
(14. 13j) and (14.161) are the results of such generalizations. In the CP-conserving limit, 
i?i 3 = i?23 = -R31 = -R32 = 0, and these expressions simplify accordingly. 

4.1 Constraints on the charged-Higgs sector 

There are three important indirect constraints on the charged-Higgs sector: the B^-B^ 
oscillations, Rf, and b — > 57. 

The mass splitting in the neutral B^ mesons is sensitive to contributions from box dia- 
grams with top quark and charged Higgs exchange [34-37], involving the Yukawa couplings 
( 13.22(1 . Indeed, the diagrams with one or two H exchanges give contributions proportional 
to (writ cot f3) 2 or (m t cot/3) 4 multiplied by functions of M#± that for large Mjj± behave 
like 1/M^ ± . These contributions to AM^ will constrain low values of M H ±, in particular 
at low values of tan/3. 

While AMs d is known experimentally to considerable precision, AM^ = (3.304 ± 
0.046) x 10~ 10 MeV [13], its theoretical understanding is more limited. The largest theo- 
retical uncertainty is related to the parameter combination j\Bs of the hadronic matrix 
element. This is only known to a precision of 10-15%. Thus, we cannot exclude models 
which give predictions for AM^ that deviate from the SM value by this order of magnitude, 
even if this deviation is large compared to the experimental precision. 

In table [T] we show the contribution to AM 6 ° that are due to the additional 2HDM 
fields, for the two parameter sets considered. It is clear that tan/3 = 0.5 is incompatible 
with the experimental and theoretical constraints on AM b °, whereas the tan/5 = 1.0 case 
is marginal. 
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tan (3 


0.5 


1.0 


2.0 


Set B: M H ± = 300 GeV 
Set A: Mjj± = 500 GeV 


270% 
140% 


45% 
24% 


10% 

5% 



Table 1: 2HDM contribution to AM 6 ° for parameter sets A and B, relative to its absolute 
value. 

As mentioned above, the b — > 57 constraints [14] also force M#± to be high, in particular 
for low tan/3. A recent analysis arrived at the bound Mjj± > 300 GeV [38]. However, at 
the very low values of tan (3 considered here, they are less severe than the AM" and R b 
constraints. The experimental constraints on Rj, (see sect. 14.51) depend on the charged 
Higgs mass as well as on the neutral Higgs spectrum. However, this constraint is for low 
values of tan (3 practically independent of the neutral spectrum. 

4.2 Higgs non-discovery at LEP 

One might think that both parameter Set A and Set B would be in conflict with the 
negative direct searches at LEP, because of the low values of M%. However, these bounds 
are marginally evaded by two facts which both dilute the experimental sensitivity. First, 
the H\ZZ coupling is suppressed by the square of the Higgs- vector- vector coupling, which 
relative to the Standard-Model coupling is 

H S ZZ : {cos (3R n + smpR j2 ], for j = 1. (4.1) 

For large values of | sinai| (which is of interest in order to maximise 7^ of (13. 24p ). Ru 

will be rather small, and the second term in (14.11) . proportional to Ru, takes over. But 

this is suppressed by the factor sin/3. For some quantitative studies of this suppression, 

which can easily be by a factor of 2 or more, see Fig. 8 in [11]. Secondly, the typical decay 

channel, Hi — > bb, is suppressed by the square of the Yukawa coupling, Eq. (13.211) . For 

small values of tan /?, this is approximately cos 2 at\ cos 2 a 2 + sin 2 (3 sin 2 a 2 . In the limits of 

interest, both terms are small. 

In the analysis of LEP data by DELPHl|], a channel-specific dilution factor C 2 is defined 

3 There are such studies also by the other LEP collaborations (see, for example [17]), but limited to the 
CP-conserving case. 
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by [16] 

&Z(h-*X) = & c zh X C Z(h^X)> ( 4 - 2 ) 

where a™ h is the Standard Model cross section for a particular Higgs mass M. In the 
2HDM, for an Hi decaying to bb or rf , the dilution is caused by the two effects discussed 
above: There is a reduced coupling to the Z boson [see fl4.ll) ] and a modified (typically 
reduced) coupling to the bb (or rf) [see (13.21 j) ] . Thus, we take 

CliH^bi) = [cos/3i2n + sin/?i? 12 ] 2 [R 2 n + sin 2 (3R\ Z \ , (4.3) 

COS l) 

and consider as excluded parameter sets those where this quantity exceeds the LEP bounds, 
roughly approximated as [16] 

C 2 z{Hi _^ bl) = 0.2 at 100 GeV, and 0.1 at 80 GeV. (4.4) 

The last term in (13.211) . involving i? 13 , is absent in the CP-conserving case. However, at 
small tan/3, it has little effect. Actually, similar results are obtained for both the bb and 
rf channels. Presumably, when these are combined, a more strict limit would be obtained. 

It is instructive to consider this expression ( 14. 3 p and the corresponding constraints in 
three simple limits: 

tan/5 <C 1 : ^ziH^bb) — cos 4 aiCos 4 a 2 <C 1. (4.5) 
This requires either a\ — > ±7r/2 or a 2 — > ±7r/2. 

tan (3 = 1: ^(Hx-^bl) — [( COSQ; i + sinai) cos a 2 ] 2 [cos 2 ai cos 2 a 2 + \ sin 2 a 2 ] <C 1. 

(4.6) 

This requires either a 2 — > ±7r/2 or ot\ — > — 7r/4 or {ai — > ±7r/2 and a 2 — > 0}. 

tan /3 3> 1 : Cz(_H"i^66) — ^ an2 ^ s * n2 ftl cos2 a 2[cos 2 cui cos 2 a 2 + sin 2 /3 sin 2 a 2 ] <C 1. 

(4.7) 

This requires <x\ — > or a 2 — > ±7r/2 or {ai — > ±7r/2 and a 2 — ► 0}. 

Furthermore, we note that at negative a\, the LEP bound is to some extent evaded for 
small and medium |a 2 | by cancellation among the two terms in the H\ZZ coupling. 
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4.3 The />-parameter constraint 

A very important constraint coming from electroweak precision data, is the precise deter- 
mination of the p-parameter [39,40]. The quantity 

A P = TFT l A ww{q 2 = 0) - cos 2 9 w Azz(q 2 = 0)] (4.8) 
M w 

measures how much the W and Z self-energies can deviate from the Standard-Model value, 
being zero at the tree level. The experiments (mostly at LEP) have put severe constraints 
on p [41]: 

p = 1.0050 ±0.0010. (4.9) 

The measured deviation from unity is accommodated within the Standard Model, and 
mostly due to the heavy top quark. 

In the 2HDM additional contributions arise [15], which are determined by the couplings 
to the W and the Z of the Higgs particles, and by the mass splittings within the Higgs 
sector, as well as the mass splittings with respect to the W and Z bosons. The simplified 
forms provided in [10] can easily be re-expressed in terms of the mass eigenvalues and 
the elements Rjk of the rotation matrix for the CP-violating basis. For the Higgs-Higgs 
contribution, we find (the relevant couplings are given in Appendix B): 

^(0)-cos 2 ^f#(0) 

= 64^2 E [{N n ^i - cos/^ 2 ] 2 + R%}F Ap (M 2 H± ,Mf) 

j 

- ^[(sin/3% - cos /3R j2 )Rk3 ~ (sm f3R kl - cos /3R k2 )R jS ] 2 F Ap {Mf, M 2 )] , (4.10) 
k>j 



where 



2 2 2 
17111712 m\ 



F Ap (ml,ml) = \{m\ + ml) ^^log^. (4.11) 

For the Higgs-ghost contribution, we have to subtract the contribution from a Standard- 
Model Higgs of mass M , since this is already taken into account in the fits, and find: 



<^(0)-cos 2 ^Aff(0) 
2 

= ^[E[ C0S ^ 1 + Sin ^ 2 ] 2 (3i 71 Ap(Mf,M 2 )-3F Ap (M^,M 2 )) 

j 
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+ 3F Ap (M^,M 2 ) -3F Ap (Mf,M 2 )j. (4.12) 

From the electroweak fits, we take M = 129 GeV, but note that this value is not very 
precise [41]. In the CP-conserving limit, these expressions (I4.10p and (I4.12p simplify con- 
siderably, since terms with i?i 3 , i? 2 3, -R31, and -R32 are absent. 

In order to keep these additional contributions (14. 1011 and (14. 12ft small, the charged 
Higgs boson should not be coupled too strongly to the W if its mass is far from those of its 
neutral partners. As a measure of the "tolerance" we take 3a, i.e., we impose \Ap\ < 0.003. 

4.4 The muon anomalous magnetic moment 

The dominant contribution of the Higgs fields to the muon anomalous magnetic moment, 
is according to refs. [42] and [25] due to the two- loop Barr-Zee effect [43], with a photon 
and a Higgs field connected to a heavy fermion loop. The contributions are given by [25] in 
terms of scalar and pseudoscalar Yukawa couplings. Re-expressed in terms of the Yukawa 
couplings of ( 13.2ip . assuming that the muon couples to the Higgs fields like a down quark, 
i.e., to $1, we find for the top quark contribution: 



(4.13) 



with N c = 3 the number of colours associated with the fermion loop, a e .m. the electro- 
magnetic finestructure constant, Q t = 2/3 and m t the top quark charge and mass, and 
m p the muon mass. The functions / and g are given in [43]. It is worth noting that the 
tan/3 factor associated with the pseudoscalar Yukawa coupling of the muon is cancelled 
by an opposite factor associated with the top quark. While the first term gives a positive 
contribution, the second one may have either sign. 

The contribution of the b quark can be obtained from ( 14. 13ft by trivial substitutions for 
Q t and m t accompanied by 

R%^Un 2 ^R%, and -L^ # - -^—R 2 jX (4.14) 

J J cos p sin p cos^p J 

in the square bracket. 

Earlier studies (see, for example, [25,42]) have focused on the contributions from rather 
light pseudoscalars and large tan/3, where the b and r contributions are enhanced by the 
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substitutions (j4.14p . At high tan /3, the 6-quark loop will indeed dominate. For small values 
of tan (3, as are considered here, the 6-quark contribution is completely negligible. 

The experimental situation is somewhat unclear, depending on how the hadronic cor- 
rections to the running of evaluated. The deviation from the Standard Model can 
be summarized as [44] 



(221 -245) x 10" 11 , e+e" 
(62 - 142) x KT 11 , T+r-, {e+e- & r+r"}, 



and represent 0.7 to 2.8 standard deviations with respect to the data [45]. Two distinct at- 
titudes are here possible. One may either fit this (positive) deviation with some new physics 
effect [25], or one may restrict new physics contributions not to exceed this contribution 
(|4.15p . We shall here follow this latter approach, and require the 2HDM contribution 
(14. 13[) to be less than 3a, i.e., |Aa M | < 300 x 10~ n . For the parameters considered here, 
tan/3 < 0{2), the 2HDM contribution to Aa^ is at most (a few) x 10~ n and therefore 
plays no role in constraining the model. 

4.5 R b 

The one-loop contributions to the Zbb coupling influence the relative branching ratio of 
Z — > bb, given by Rb, which is known to 0.05%, or 1.25 MeV precision [13]. In the SM there 
are significant contributions proportional to mf . In the 2HDM there are additional one- loop 
contributions due to triangle diagrams involving charged and (non-standard) neutral-Higgs 
fields. For the CP-conserving case, these were given in [46]. In the general CP- violating 
case, the charged-Higgs contribution, Eq. (4.2) of [46], remains unchanged, but we find 
that the neutral-Higgs part, Eq. (4.4), gets modified to 

967rsm 9 w cos 2 9 w 77% 



k=l 



tan Q 

- (sin[3R kl - cos /3R k2 )R j3 ] £ R n R k3 p 4 (m 2 z , M 2 , M 2 , 0) 

cos p 

R ■ 

- (cos f3Rj! + sin (5R j2 )^^ pi(m 2 z , M 2 , m|, 0) 



cos/3 " 1 
f R 2 

2Q b sin 2 9 W (1 + 2Q b sin 2 9 W ) ( + tan 2 (3R% ) p 3 (m 2 z , M 2 , 0) 

COS l) 
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+ 2Q b sin 2 9 W (1 + 2Q b sin 2 6 W ) p 3 (m 2 z , m|, 0) J. (4.16) 

The functions p 3 and P4 are various combinations of three-point and two-point loop in- 
tegrals [46]. For the numerical studies, we use the LoopTools package [47,48]. Again, 
this expression (I4.16P is more complicated than those of the CP-conserving limit, but the 
additional terms have little quantitative importance at low tan (3. 



tan j3 


0.5 


1.0 


2.0 


Set B: M H ± = 300 GeV 
Set A: M H ± = 500 GeV 


5.6 
3.3 


1.4 
0.83 


0.35 
0.21 



Table 2: 2HDM contribution to R b for parameter sets A and B, in units of the uncertainty, 
a(R b ) = 1.25 MeV. 

At small tan/?, the charged-Higgs contribution, which behaves like (m 4 /tan/3) 2 due 
to the H + bi and H~bt couplings, dominates. For M#± = 500 GeV (Set A) and Mh± = 
300 GeV (Set B) the 2HDM contributions to R b are given in Table [21 For the two larger 
values of tan/3, this is compatible with the experimental uncertainty, whereas for tan/? = 
0.5 it amounts to a substantial conflict, especially for the lower value of M#± (Set B). The 
neutral-Higgs contribution, as given by (14. 16ft . is smaller, by three orders of magnitude. 

5 Overview over allowed parameters 

5.1 Variations of mass parameters around Set A 

We start this discussion of model parameters by a survey of how the allowed regions of the 
«2~«3 space depend on the mass parameters, in particular M3, M#± and p 2 . It turns out 
that while stability is readily satisfied for "relevant" mass parameters, unitarity excludes 
sizable regions of parameter space. Ignoring experimental constraints, a low-mass spectrum 
is in general easier to accommodate than one where some Higgs particles are heavy. In 
many cases, non-zero values of \$ and A7 also have a tendency to reduce the allowed 
parameter space. 

In Fig. [T]we show for Set A the allowed regions in the a2-«3 plane, for a few represen- 
tative values of a±, focusing on stability (or positivity) and unitarity. For the considered 
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Figure 1: Allowed and forbidden regions in the «2^«3 plane, for tan/3 = 0.5, 1 and 2, 
parameter Set A and selected values of a\. Green (G): stability, unitarity, direct search 
and Ap constraints satisfied; red (R) [purple (P)]: stability, unitarity satisfied, but direct 
search [and Ap] constraints not satisfied; yellow (Y): stability (but not unitarity) satisfied; 
white: stability violated. 



parameters, much of the 02-0(3 plane is actually in violation of stability, as shown in white. 
Next, there are significant areas (yellow) where stability is satisfied, but (tree-level) uni- 
tarity is violated. Finally, in darker colours, we indicate where also unitarity holds. Part 
of these areas (typically, small \a 2 \, as discussed in sect. 14. 2p are in conflict with the direct 
search (LEP) data, and indicated in red. Some areas (purple, typically, larger |a 2 |) violate 
the Ap constraint, whereas the remaining areas in green give viable models. The symmetry 
given by Eq. (13.111) is evident in the figure: the case of tan (3 = 2 can be obtained from the 
case of tan (3 = 0.5 by suitable reflections, apart from the experimental constraints, which 
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in the 2HDM (II) distinguish tan/3 and cot/3. For tan/3 > 1, the direct search constraints 
shift from those of ( 14.5ft towards those of ( 14.71) . 

While this figure can not be directly compared with Fig. 7 of ref. [11], since we here 
keep M3 fixed, it was found that the unitarity constraints of [30,32] are more restrictive 
than the order-of-magnitude estimate adopted in [11]. It should also be noted that with 
Im A5 = 0, as given by the parameter Set A, Im X G and Im A7 will be non-zero in the general 
CP-violating case. 
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Figure 2: Physically allowed regions in the 02-03 plane, for tan/3 = 0.5 and 1, and selected 
values of a%. Left: Variations of M3 around parameter Set A. Blue: M3 = 500 GeV (Set A), 
vertical lines: M 3 = 550 GeV, yellow: M 3 = 600 GeV. Note that some regions overlap. 
In particular, green denotes regions allowed for both M 3 = 500 GeV and M 3 = 600 GeV. 
Right: Variations of M H ± around parameter Set A. Yellow: M H ± = 300 GeV, vertical 
lines: M H ± = 400 GeV, blue: M H ± = 500 GeV (Set A). 
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Dependence on M 3 and M 2 . The dependence of the allowed regions on M 3 , the heaviest 
neutral Higgs boson, is illustrated in Fig. [2] (left), for the other parameters kept fixed at 
the values of parameter Set A. The figure shows the allowed regions for M 3 = 500 GeV 
(blue, default of Fig. [T]), 550 GeV (vertical lines) and 600 GeV (yellow). Smaller allowed 
regions are also found at 450 and 650 GeV (not shown), but nothing neither at 400 GeV 
nor at 700 GeV. 

As M2 approaches Mi, there are still regions where stability is satisfied, but unitarity 
is only satisfied in very small regions. Similarly, as discussed above, when M 2 approaches 
M 3 , the allowed regions tend to be restricted to small values of |a 2 | - > 0. Thus, by (13.341) . 
there is in this limit of Mi -C M 2 ^ M 3 , little CP violation in the top-quark sector. 

Dependence on M^. The dependence of the allowed regions on Mh± , the charged Higgs 
boson, is illustrated in Fig. [2] (right), for the other parameters kept fixed at the values of 
parameter Set A. The allowed region is seen to increase a bit for lower values of M H ±, 
but such values are constrained by the AM^, Rb [13] and b — > 57 data [14], especially 
at low values of tan f3. On the other hand, the allowed regions shrink for high values of 
Mff± . Nothing is allowed at > 600 GeV. The shrinking of the allowed regions for high 
values of M H ± is due to the Ap constraint at large |ck 2 | and/or negative ati, as well as the 
direct search (LEP) constraint at small |ct 2 1 • Eventually (like for high values of M 3 ) also 
the unitarity constraint excludes high values of M^. 

Dependence on fi 2 . Since A4 + Re A5 is bounded from below by the stability requirement 
[see (IA.23I) ]. the parameter /1 2 normally plays a role in pushing M H ± to values that are 
high enough to evade the i4, AMB d [13] and b — > sj constraints. In Fig. [3] we show how 
the allowed region in the 02-0:3 plane shrinks for "low" and "high" values of fi 2 , when 
the Higgs boson masses and other parameters are kept fixed. For the considered spectrum 
of masses, a range of negative values of /1 2 is allowed (for fi 2 = —(400 GeV) 2 there are 
no allowed regions). For increasing positive values of /1 2 , the allowed regions shrink away 
before fi 2 = (400 GeV) 2 (which is not allowed). Of course, these critical values depend on 
the mass spectrum adopted in parameter Set A. 
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Figure 3: Physically allowed regions in the a.2-013 plane, for tan/? = 0.5 and 1, and selected 
values of a\. Variations around parameter Set A; /i 2 < and /i 2 > shown separately. 
li 2 < 0: Yellow: li 2 = 0, green: /1 2 = -(100 GeV) 2 , vertical lines: /1 2 = -(200 GeV) 2 . 
li 2 > 0: Yellow: /x 2 = (100 GeV) 2 , blue: /i 2 = (200 GeV) 2 , vertical lines: li 2 = (300 GeV) 2 . 
Note that some regions overlap. In particular, for /1 2 > 0, green is allowed for both 
il 2 = (100 GeV) 2 and ll 2 = (200 GeV) 2 . 



5.2 Variations of mass parameters around Set B 

We shall here briefly review the light-Higgs scenario of parameter Set B, but recall that 
the tan/5 = 0.5 case is essentially ruled out by the AM 6 ° data. 

In Fig. H] we show for li 2 = how stability can be satisfied in the whole a^-ctz plane. 
For tan = 1, stability and unitarity would actually allow any value of the rotation matrix 
R, i.e., any values of a%, ot 2 and a 3 . However, except for this special case, unitarity is only 
satisfied in parts of the plane. 

The direct search (LEP) constraint severely cuts into the allowed regions of this light- 
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Figure 4: Stability and unitarity in the a 2 -a^ plane, for tan/3 = 0.5, 1 and 2, parameter 
Set B with M 3 = 400 GeV, fi 2 = and selected values of a±. Colour codes as in Fig. [H 



Higgs (set B) scenario, with C 2 = 0.1 [16], see sect. 14.21 As noted there, this constraint 
tends to exclude the central range of small \a 2 \ and positive and small at\. 

Increasing fi 2 to (200 GeV) 2 , the most striking change is perhaps the emergence of large 
regions where stability is not satisfied (indicated in white in Fig. Another interesting 
observation is that for tan/3 = 1.0 and negative values of «i, the picture is little changed 
from the case of [i 2 = 0. 

With a higher value of M3 (M3 = 600 GeV, but fi 2 = 0), the unitarity constraint 
excludes most of the «2-«3 plane. For tan/3 = 1.0 the whole plane is excluded. 



27 



tan (3 = 0.5 tan /3 = 1.0 tan /3 = 2.0 




Figure 5: Stability and unitarity in the «2-«3 plane, for tan (3 = 0.5, 1 and 2, parameter 
Set B with M 3 = 400 GeV, fi 2 = (200 GeV) 2 and selected values of a%. Colour codes as in 
Fig.ffl 

5.3 Non-zero values of ImAs, ReA6, ReXj 

This subsection will present a brief discussion of how the allowed regions get modified 
for non-zero values of parameters which are normally set to zero, in order to control the 
amount of Flavour-Changing Neutral Currents. Thus, adopting a non-zero value for any 
of them in a "realistic" model would have to be done with an eye to these effects. 

5.3.1 Non-zero values of ImA 5 

The cases of positive and negative values of Im A5 can be related. This is seen as follows: 
by the transformation C3 of (13.61) and (13.81) . the mass-squared elements J\A 2 3 and .M23 
flip signs, and the auxiliary quantities ImAg and IH1A7 change signs without altering the 
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stability or unitarity constraints. The only effect will be that CP-violating effects change 
sign. Thus, we may restrict the discussion of non-zero ImAs to ImAs > 0. However, the 
full range (13.7p of 03 now has to be considered. 

We show in Fig. [6], for some values of tan/3 and a±, and for parameter Set A how 
non-zero values of ImAs also provide allowed regions in the 02-03 plane. But these are 
typically smaller for higher values of ImAs, an d rather scattered. However, they can lead 
to significant CP violation, since allowed ranges of \a2\ tend to be at intermediate values 
of \a 2 \ [see Eq. (15T71V . 

Varying lm A,-, 
tan 5=0.5 tan/? = 1.0 tan = 2.0 



jt/2 
an 



n/2 



■jt/2 



■jt/2 



■jt/2 



?r/2 



-tt/2 





4 






2' 











-¥72 



ttJ2 -i/2 



a 


1 




2 














f 












c 


1 


































*} 




















U 








a 


: = 































a 


1 = 
>| 


B 


r/6 
r 






m 












a 


1 = 




r/3 




4 


r 

































_ 




"I 


1- 
























































































































































































■s/2 „ 


jt/2 



Figure 6: Physically allowed regions in the 02-^3 plane, for tan/3 = 0.5, 1 and 2, and 
selected values of o.\. Variations of ImAs around parameter Set A. Blue: ImAs = 0, 
vertical lines: Im A5 = 2, yellow: Im A5 = 4. 
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5.3.2 Non-zero values of ReA fi or ReA 7 



Up to this point, we have kept ReAg = ReA7 = 0. However, we have treated ImAg and 
IH1A7 as auxiliary quantities derived from the spectrum, the rotation matrix and ImAs 
via Eq. f 1 3 . 3 H . In general, they will be non-zero. This might lead to too large flavour- 
changing neutral couplings, due to the violation of the Z 2 symmetry [2,27,28]. We have 
not investigated this constraint quantitatively, but note that in many cases the imaginary 
parts of Ag and A7 can be shifted to the imaginary part of A5 [see ( 12. 10p ] . This will however 
lead to a modification of the rotation matrix R and/or the spectrum, by for example a 

shift in M3 according to the approach of [11]. 

tan/9 = 0.5 tan /? = 1.0 tan ft = 2.0 
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Figure 7: Physically allowed regions in the «2-«3 plane, for tan/? = 0.5, 1 and 2, and 
selected values of ol\. Parameters correspond to Set A, except that ReA 6 = —1. Colour 
codes as in Fig. [TJ 
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The Yukawa interactions (see sect. 13. 5j) couple the Higgs fields to a left-handed doublet 
and a right-handed singlet quark field. However, these need not be in the flavour basis in 
which the mass matrices are diagonal. The Z2 symmetry, which is imposed to stabilize 
Model II [2,28] is broken by the and Im A5 terms, as well as by the X% and A7 terms. 
However, one may adopt the attitude that these terms, which arise naturally in the MSSM 
[18-22] are constrained and subdominant. Additionally, one may argue that the FCNC's 
are suppressed by powers of the quark masses [49,50], thus evading the experimental 
constraints involving the first two fermion generations. 

We shall here only discuss the case of either Re Xq or Re A7 being non-zero, the other 
being zero. Because of the symmetry (13.111) . (13.1 2p . we may restrict this discussion to 
ReA 6 7^ 0, Re A 7 = 0. Analogous results for ReA 6 = 0, ReA 7 7^ can be obtained from 
these by the inversion tan/3 <-> 1/ tan/? according to (13.111) and (13.121) . 

Even though Re Xq (or Re A7) significantly different from zero may turn out to be ruled 
out by the constraints on FCNC's, we find it instructive to see how the otherwise allowed 
regions change when these parameters are introduced. 

In Fig. [7] we show the "allowed" regions corresponding to parameter Set A, except 
that ReAe = —1. The allowed regions are qualitatively rather similar to those of Fig. [U 
However, with ReA6 = +1, stability tends to be violated in most of the 0:2-03 plane. Also, 
when ReA 6 decreases to —2 or —3, there is nothing allowed at tan/? = 1 and tan/3 = 0.5, 
respectively. 

6 CP violation in tt production at the LHC 

In order to illustrate the CP-violating effects that can be observed at the LHC, resulting 
from mixing in the Higgs sector, we consider the process 

pp^ti+X, (6.1) 

which at high energies is dominated by the underlying process 

gg -> tt. (6.2) 
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In the presence of CP violation, correlations will then be induced at the parton level 
involving the t and t spins, denoted S\ and 82, and their cm. momentum p [8]: 

(p-(s 1 -s 2 )}, and (p-(«iXs 2 )). (6.3) 

These correlations are determined by the CP- violating combination (13.231) of Yukawa cou- 
plings, multiplied by certain loop integrals, and convoluted over the gluon-gluon cm. 
energy [8, 11]. 

The t and i quarks decay fast enough that hadronization effects do not smear out these 
CP-odd correlations. They can thus be treated perturbatively. Consider the semileptonic 
decays: 

t -> l + vib, i -> ruib, (6.4) 

with I = e~ or jjT . The lepton energies will inherit an asymmetry from the first correlation 
given in Eq. (16. 3p . accessible via the observable: 

A 1 = E+-E_. (6.5) 

An important question is whether or not this can be large enough to be measurable. 

In order to have a significant observation, the expectation value (Ai) must compare 
favourably with the statistical fluctuations, which behave like y/N, where N is the number 
of events. In order to assess this, it is convenient to consider the "signal to noise" ratio [8] , 

I = , W ■ (6.6) 

where the denominator gives the statistical width of the observable (16.51) . 

Other limitations of this approach include the need for very good lepton energy calibra- 
tion, and the assumption of no (anomalous) right-handed couplings in the thW coupling 
(see [51]). 

We display the ratio S/N in Fig. [HI as a function of the mass of the lightest Higgs 
boson, Mi, for two rotation matrices R, given by the angles {«i, a^}, and for two 
values of tan/3. These parameters are chosen such that the model is consistent (stability 
and unitarity) and satisfies the experimental constraints (see the discussion in Sect. 15.11 
and Fig. [T]). The actual evaluation of this quantity (16.61) follows the approach of [11], using 
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the LoopTools package for the loop evaluations [47,48] and the CTEQ6 parton distribution 
functions [52] to describe the gluon content of the proton. 
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Figure 8: Signal to noise ratio vs. Mi, otherwise parameter Set A. Rotation matrix R 
defined by angles {ai,«2,«3} as indicated. Actually, in the left panel, we show S/N for 
—A\. Dotted curves: contribution of the lightest Higgs boson only. 

While the allowed regions, for tan/3 = 0.5, in Fig.[T]only are rather small and scattered, 
they remain allowed as Mi is increased towards M 2 . For these two rotation matrices, the 
contribution of the lightest Higgs boson, Hi, is actually the same (apart from the over-all 
sign), as shown by the dotted curves. In terms of the quantities (13.231) . \jq P \ = 1.936 is 
the same for both cases. However, for the case shown on the left, the contribution of if 2 
(with M 2 = 300 GeV), is with ^opHcp = -°- 75 

more efficient in reducing the effect of 
Hi than for the case on the right, where the corresponding ratio is —0.44. 

In general, the values fall with increasing Mi, since the loop integrals decrease. However, 
there is a resonance in one diagram at Mj = 2m t ~ 350 GeV (for more details, see [11]), 
this is the reason for the increase beyond 250 GeV. 

In Fig. [9] we show, for the same rotation matrices as in Fig. [Hi how S/N varies with 
M 2 . Here, some of the curves are cut off at low or high values of M 2 , since the model 
becomes inconsistent or experimentally excluded, as discussed in Sect. 15.11 In some cases, 
an enhanced negative interference occurs as M 2 ~ 2m t . 

In this study, the parameter \j? was allowed to float, as compared with parameter Set A, 
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Figure 9: Signal to noise ratio vs. M 2 , otherwise parameter Set A. Rotation matrices R as 
in Fig. El Dashed line: value of M 2 in parameter Set A. 

in order to extend the range of allowed values of M 2 . In all cases, except the right panel, 
with tan/3 = 1, the value of p 2 had to be adjusted as follows: At low values of M 2 , a lower 
value of jj 2 was needed, and at high values of M 2 , a higher value was needed. 

Finally, in Fig. [TUJ we show the variation of S/N with M 3 , keeping Mi and M 2 fixed. 
There is a tendency for the value to increase with M3 (because of reduced destructive 
interference), but variations of M3 are only allowed within a rather restricted range (see 
again Sect. 15.11 and Fig. [2]). 



7 Summary and conclusions 

The constraints of stability and tree-level unitarity exclude most of the multidimensional 
2HDM (II) parameter space. Furthermore, the direct searches (LEP) [16] and the Ap [13] 
constraints exclude certains domains of the parameters. Finally, the direct searches, as well 
as the AM Bd , R b [13] and b — > 57 [14] constraints exclude a light charged Higgs boson, in 
particular at low values of tan j3. The remaining pockets of allowed parameters will for low 
values of tan j3 allow CP violation in the top-quark sector, due to the exchange of Higgs 
bosons that are mixed with respect to CP. 

In order to maximize the CP violation in the top-quark sector that might be measurable 
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Figure 10: Signal to noise ratio vs. M 3 , otherwise parameter Set A. Rotation angles 
{ai, CK2, 0(3} as indicated. Dashed line: value of M 3 in parameter Set A. 

at the LHC, we focus on parameters where (i) the lightest Higgs boson is rather light, in 
order to maximize the relevant loop integrals, and (ii) where the product of the CP- violating 
Yukawa couplings, parametrized by [see Eq. ( I3.23p ] are large. The latter constraint 
requires tan f3 to be small, and | sin a,\ sin 2a 2 1 to be large. 

In summary, we note that even in the face of a variety of experimental constraints, the 
model is consistent in a number of regions in parameter space. Apart from exceptional 
points, these allowed regions yield CP violation in ti final states produced at the LHC, at 
a level which can be explored with a data sample of the order of 10 6 semileptonic events. 
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Appendix A 

In this Appendix we describe the method used for formulating the necessary and sufficient 
conditions for the potential to satisfy stability. For earlier approaches to stability, all 
restricted to simpler potentials, see [53-56]. We start by rewriting the Higgs doublets as: 



$! = ||$i||$i, 



$5 



|$2||$2 



(Al) 
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where is the norm of the spinor $j, and $j is a unit spinor. By SU(2) invariance, 

only four combinations of fields may appear: 

$t $1 = 11$^ $ t $2 = ||$ 2 ||2 j = || $l || . || $2 || fat $ t $2 = [$t $l ]* 

We let the norms ||$j|| of eq. ( 1A.1I) be parametrized as follows: 

| |$i| | = r COS7, ||$2|| = TSU17. (A. 2) 

The complex product $2 ' ^1 between the two unit spinors will be a complex number x + iy 
with \x + iy\ < 1. 

Using this parametrization, we can write: 

= r 2 cos 2 7, $2$2 = r 2 sin 2 7, 
$2^1 = r 2 cos 7 sin 7(2 + iy), §\§ 2 = r 2 cos 7 sin 7(2 — iy), (A. 3) 

where r > 0, 7 G [0, 7r/2] and x 2 + y 2 < 1. 
The potential can now be written as 

V = r 4 V 4 + r 2 V 2} (A.4) 

with the quartic part: 

V4 = Ai At + X 2 A 2 + A3A3 + A4A4 + Re A 5 A 5 + Im X 5 A 6 

+Re A 6 A 7 + Im \ 6 A 8 + Re A7A9 + Im X 7 A 10 (A.5) 

where 

Ai = ^cos 4 7, A 2 = isin 4 7, A 3 = cos 2 7 sin 2 7, 

At = (2 + y ) cos 7 sin 7, A 5 = (x — y ) cos 7 sin 7, A 6 = 2xy cos 7 sin 7, 

= 2x cos 3 7 sin 7, A 8 = 2y cos 3 7 sin 7, 
v4 9 = 2x cos 7 sin 3 7, A 10 = 2?/ cos 7 sin 3 7. (A. 6) 

By introducing polar coordinates x = pcos9 and y = psin.9, these coefficients can be 
written 

A x — - cos 4 7, A 2 = - sin 4 7, A 3 = cos 2 7 sin 2 7, 
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At = p 2 cos 2 7 sin 2 7, A 5 = p 2 cos 2 7 sin 2 7 cos(20), A Q = p 2 cos 2 7 sin 2 7 sin(26>), 
A 7 = 2p cos 3 7 sin 7 cos 9, A s = 2p cos 3 7 sin 7 sin 9, 

A 9 = 2p cos 7 sin 3 7 cos 0, A w = 2p cos 7 sin 3 7 sin 0, (A. 7) 

where 7 G [0, vr/2], p G [0, 1] and 9 G [0, 2tt). 
A.l Symmetries of V4 

We note some symmetries of the quartic potential under the parametrization (1A.2jl . They 
are conditional symmetries, where the potential is invariant with some "compensating" 
interchange of As: 
Symmetry I: 
The interchange 

y <-> -y or 9 ^ 2n - 9 (A.8) 

can be compensated by 

{IrnAs, ImA6, Im A7} <-> {— Im A5, — ImA6, —lm\ 7 } (A. 9) 

This is of course nothing but the reality condition of the potential. 

Symmetry II: 

Under 

7^ I -7 (A.10) 

together with 

Ai <-> A2, and A6 <-> A7 (A. 11) 

the potential is invariant. This is the symmetry under the interchange $1 <-> $2- 

Symmetry III: 

Under 

(x, y) <-> (— — y) or <-> tt + 9 mod 2tt (A. 12) 

together with 

{A 6 ,A 7 }^{-A 6 ,-A 7 } (A.13) 
the potential is invariant. This is related to the well-known Z 2 symmetry [2,28]. 
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Symmetry IV: 

Under 

(x, y) <-» (y, x) or 9^> — — 9 mod 27r (A. 14) 

together with 

ReA 5 ^ — ReA 5 and {Re A 6 , Re A 7 } <-> {Im A 6 , Im A 7 } (A. 15) 

the potential is invariant. 

A.2 Stability 

For the stability condition to be satisfied, V4 must be positive for all combinations of 
7 G [0, 7r/2], p G [0, 1] and 9 G [0, 2rr). This is both a necessary and a sufficient condition. 
Whenever A4 < 0, the potential will have its global maximum and minimum when p = 1. 
Thus, it is sufficient to check that ^(7, 9; p = 1) satisfies stability when A4 is non-positive. 
In order to see this, we return to eqs. (1A.5I) and ( 1A.6I) and rewrite the potential as 

V4 = A 4 (x 2 + y 2 ) cos 2 7 sin 2 7 + y) 

where h(x, y) is a harmonic function. The maximum principle tells us that h(x, y) will 
attain its global minimum at the boundary where x 2 + y 2 = 1 (p = 1). Whenever A4 < 0, 
the term A 4 (x 2 + y 2 ) cos 2 7 sin 2 7 will also attain its minimum whenever p = 1, and so will 

Some points from the parameter space give us some rather simple stability conditions. 
We now turn our attention towards these special points. 
7 = or 7 = 7r/2 

First we consider the boundary points 7 = and 7 = tt/2. 

^4(7 = 0) = ^ 
14(7 = tt/2) = y 
This leaves us with the "trivial" stability conditions of (13.201) : 

Ai > and A 2 > 0. (A.16) 
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p = 

Considering the points where p = 0, we find that 

Vi(p = 0) = — - cos 4 7 + sin 4 7 + A3 cos 2 7 sin 2 7. (A. 17) 

Thus, we must require that 

A 3 > ~\ ( r4- + A 2 tan 2 7 ) . (A.18) 
2 \tan 7 / 

The right-hand side has its minimum for tan 2 7 = A/A1/A2, and we obtain the following 
necessary constraint on A3: 

A 3 > -\A^, (A.19) 

in agreement with [53]. 

A.3 The limit A 6 = A 7 = 

It is instructive to consider the simple limit 

A 6 = A 7 = (A.20) 
Then the quartic part of the potential can be written 

V A =y cos 4 7 + y sin 4 7 + [A 3 + p 2 (A 4 + Re A 5 cos 29 + Im A 5 sin 29)} cos 2 7 sin 2 7 (A. 21) 
This expression has the same structure as (1A.171) . with 

A3 -> A 3 + p 2 [A 4 + Re A 5 cos 29 + Im A 5 sin 29] . (A.22) 
Thus, the condition for stability can be adapted from ( 1A.19I) . leading to: 

A 3 + min[0,A 4 - |A 5 |] > -v 7 ^, (A.23) 
as obtained by [53]. However, we stress that this constraint only applies when X e = X 7 = 0. 
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Appendix B 



The Higgs-vector-boson couplings can be extracted from the covariant derivatives: 
(W 1 )t(L»^ 1 ) + (D^ 2 )\D^ 2 ) 

-Z^{- cos 29 w (ipid IJi <pT - vid^Vi + ^2 d ^t ~ vt d ^ 2 ) 



2 cos 9]y 

+ ^^"[fai + iXi)d^T - tPid»(Vi + + (^2 + %)<9 M v?2 - ^2 ^(^2 + %)] + 



(B.l) 



These Higgs fields have to be transformed into the physical basis: 

■ ni = R ji H j , 1 = 1,2,3, (B.2) 

with 



(B-3) 



(pf\ _ /cos 13 -sin/A /G ± \ /^A _ /cos/3 -sin /ft /G° 
~ V sin/3 cos/3 J ' VW ~ U 11 / 9 c 08 / 9 / \V3 

With all momenta incoming (in an obvious notation), we find 

7TT+U- cos 20 w + 

2 COS V\y 

-ig cos28 w . 
ZG G • 2cos0 w ^-Pm)' 

Zff^fc : ~ g , [(sin/3% - cos PR j2 )R k3 - (sinf3R kl - cos /3R k2 )R j3 M ~ it), 

Z COS (7^ 

: ^hr ( cos ^ 1 + sin ^)(pi - P j), (b.4) 

z cos t>w 

and 

W^H^Hj : |[±i(sin/3% - cos/3%) + %](p> - P J), 
W ± G^H j : ^(cos/3% + sin/3%)(p£ -pj), 

iy ± G =F G' : f(pP-p?) (B.5) 

There are no ZH ± G T or W ± H Zf G° couplings. The CP-conserving limit is obtained by 
evaluating R for a 2 = 0, a 3 = 0, a\ = a + n/2, with the mapping H 1 ^ h, H 2 ^ —H and 
H 3 — > A. In that limit, we recover the results of [10]. 
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